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We derive the Hamiltonian in the rotating frame for NMR quantum computing with homonu- 
cleus molècules as its computational resource. The Hamiltonian thus obtained is different from 
conventional Hamiltonians that appear in literature. It is shown that control pulses designed for 
heteronucleus spins can be translated to pulses for homonucleus spins by simply replacing hard 
pulses by soft pulses with properly chosen pulse width. To demonstrate the validity of our Hamil- 
tonian, we conduct several experiments employing cytosine as a homonucleus molecule. All the 
experimental results indicate that our Hamiltonian accurately describes the dynamics of the spins 
and that the conventional Hamiltonian fails. Finally we use our Hamiltonian for precise control of 
field inhomogeneity compensation with a pair of 7r-pulses. 
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I. INTRODUCTION 



Quantum computation currently attracts a lot of at- 
tention since it is expected to solve some of computa- 
tionally hard problems for a conventional digital com- 
puter |]| . Numerous realizations of a quantum computer 
have been proposed to date. Among othcrs, a liquid-state 
NMR (nuclear magnètic resonance) quantum computer is 
regarded as most successful. Early experiments demon- 
strated quantum teleportation 0], quantum search algo- 
rithm Q , quantum error correction Q , and simulation of 
a quantum mechanical system p| . Undoubtedly, demon- 
stration of Shor's factorization algorithm 6] is one of the 
most remarkable achievements in NMR quantum compu- 
tation. Although the number of admissible qubits in a 
liquid-state NMR quantum computer is suspected to be 
limited up to about ten due to poor spin polarization at 
a room temperature, a liquid-state NMR quantum com- 
puter is one of few quantum computers that are capable 
of running nontrivial quantum algorithms thanks to well 
established NMR technology. 

Since the number of qubits within heteronucleus spins 
is practically limited to two or three, the use of homonu- 
cleus spins is inevitable if we try to equip an NMR with 
a large number of qubits. It should be pointed out, how- 
ever, that liquid-state NMR of homonucleus molècules is 
still poorly understood and literature dealing with this 
subject often lacks sòlid ground. Although the prod- 
uct operator formalism Q has been extensively employed 
to implement quantum algorithms with an NMR quan- 
tum computer, people overlooked significance of the gen- 
uine Hamiltonian. Actually there is a subtle difference 
between the conventionally used Hamiltonian and the 
proper Hamiltonian for homonucleus spins. It is, there- 
fore, urgently required to establish theoretical foundation 
underlying a liquid-state NMR quantum computer with 



homonucleus molècules. 

Suppose we would like to implement a quantum al- 
gorithm whose unitary matrix representation is E/ a ig- If 
the Hamiltonian H depends on the control parameters, 
which we write collectively as "f(t), the time evolution 
operator is given by 



£7[ 7 (í)]=Texp 



H{l{t))dt 



(D 



where T stands for the time-ordering product. We use 
the natural unit in which K = 1. Optimal control of the 
quantum computer requires a control function j(t) that 
produces the specified quantum algorithm Í7[7(í)] = C/ a i g 
in the shortest possible time T. Recently, numerical 
scheme to find the optimal control has been worked out 
for fictitious Josephson junction qubits, where polygonal 
paths in the parameter space has been utilized [g, 0- 
For time-optimal control of an NMR quantum computer, 
anothcr method employing the Cartan decomposition of 
SU(2 n ) has been proposed [l(| and has been demon- 
strated experimentally with a two-qubit heteronucleus 
molecule [ÏJJ. We note that exact optimal control has 
been found for holonomic quantum computation in an 
idealized situation [ï^ |. 

This paper has three aims: (1) to provide the theo- 
retical foundation for an NMR quantum computer with 
homonucleus molècules, (2) to show that any pulse 
sequence designed for heteronucleus molècules can be 
translated into that for homonucleus molècules, and (3) 
to demonstrate experimentally that our Hamiltonian ac- 
curately describes the dynamics of the spins. For these 
purposes, we carefully cxaminc the Hamiltonians for 
NMR spin dynamics. Although the Hamiltonian for a 
homonucleus molecule is the same as the one for a het- 
eronucleus molecule in the laboratory frame, the former 
looks quite different from the latter in a rotating frame. 
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This paper is organized as follows. In section II we 
study Hamiltonians of homonucleus as well as heteronu- 
cleus molècules. We carefully examine how they are 
transformed in a rotating frame and what is appropriate 
approximation to be employed. Surprisingly, our result- 
ing Hamiltonian is difFerent from the convcntional Hamil- 
tonian. In section III, we conduct several experiments to 
verify our analysis by taking cytosine as an example of 
homonucleus molècules. We execute the Deutsch-Jozsa 
algorithm, execute the pulse sequence for pseudo-pure 
state preparation, and verify the robustness of two-qubit 
cntangling operations. As an application of the correct 
form of the Hamiltonian we implcmcnt field inhomogene- 
ity compensation using a pair of 7r-pulses in section IV. 
Section V is devoted to conclusions and discussion. 



II. HAMILTONIAN IN ROTATING FRAME 

In this section, we write down the Hamiltonian of spin 
dynamics in the laboratory frame and transform it to 
the one in a rotating frame. Although the Hamiltonian 
for a homonucleus molecule has the same form as the 
one for a heteronucleus molecule in the laboratory frame, 
Hamiltonians in a rotating frame differ from each other. 

We restrict ourselves within two-qubit molècules for 
simplicity. Generalization to molècules with more qubits 
is straightforward. As an example of heteronucleus 
molècules, we refer to 13 C-labeled chloroform. The qubits 
are spins of 13 C and H nuclei. We take cytosine solved in 
D2O as an example of homonucleus molecule. The qubits 
are spins of two hydrogen nuclei (protons) in this case. 



A. Heteronucleus molecule 

1. Experimental setup 

A liquid-state NMR consists of three parts as described 
in |l3j . The first part is magnètic coils; a superconduct- 
ing coil to generate a homogeneous static magnètic field 
and a normal conducting coil to generate temporally con- 
trolled field gradients. The second part contains reso- 
nance circuits for applying radio frequency (rf) magnètic 
fields to the sample. They are also used to pick up rf 
signals from the sample. The third part is an assembly 
of electrònic circuits to feed rf pulses into the resonance 
circuits and to detect the signals picked up by the coils. 

The NMR setup for heteronucleus molècules is shown 
schematically in Fig. ^ The first and second spins 
have respective resonance freqüències ujq,í (i — 1,2), 
which are also called Larmor freqüències. Their reso- 
nance freqüències are widely different for heteronucleus 
molècules under consideration. Hcnce two sets of res- 
onance circuits and assembly of electrònic circuits are 
required. The large difference of the resonance freqüèn- 
cies, Alüq = UJ0.2 — ^0,1, allows us to address each spin 
individually with a short pulse. 
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FIG. 1: NMR setup for heteronucleus molècules. 



The oscillator i (i — 1,2) in the third part generates 
an rf elèctric wave with frequency Wrfj. The sequencer 
i modulates the rf wave to shape a designed pulse. A 
typical temporal duration of a pulse, which is called the 
pulse width, is of the order of 10 fis. The rf pulses are 
amplificd and fed into the resonance coil i, which gener- 
ates rf magnètic fields applied to the sample in the test 
tube. Precession of spins in molècules appears as rota- 
tion of magnetization of the sample and induces a signal 
at the coil i. The receiver i detects the signal. The direc- 
tional coupler prevents transmission of the rf pulse from 
the amplificr to the receiver. 



2. Heteronucleus molecule in rotating frame 
The two-qubit Hamiltonian in the laboratory frame is 

H = H + H r f,i + H r ifl- (2) 
Here the system Hamiltonian Ho is defined as 0] 

Hq = — CÜQlI- 



I - lü 0i2 I <g> I z + Jh®h, (3) 

k— x,y ,z 



where Ik = Cfc/2, (Jk being the /c-th Pauli matrix, and 
I is the unit matrix of dimension two. The first two 



3 



terms in Hq describe free precession of the spins in a 
static magnètic field while the third term describes the 
intramolecule spin interaction with coupling strength J. 

H T í t i (i = 1, 2) represents the action of the rf magnètic 
field generated by the coil í and hcncc is called the control 
Hamiltonian. Thcir cxplicit forms are 

ií rfi i = -2w M cos(w rU í - (j)i)(I x ®I + gI® I x ), (4) 
H x f í2 = -2wi i2 cos(w r f i2 í - ^(g^Ix &> / + 1 <&I X ).(5) 

Here, the amplitude of the rf pulse u>i t i, the frequcncy 
of the pulse ui r f a and the phase of the pulse 4>i are con- 
trollable parameters. We may assume, without loss of 
generality, that the rf field is applied along the x-axis in 
the laboratory frame. In the above equations we intro- 
duced the ratio of resonance freqüències of the two nuclei, 

9 = —- (6) 

We shall examine the transformation law of the Hamil- 
tonians from the laboratory frame to a rotating frame. 
The spin dynamics in the laboratory frame is governed 
by the Liouville equation 

i% = [H,p], (7) 

where p is the density matrix of the system under con- 
sideration. The unitary operator 

U = exp(-íw rot ,iJ z í) ® exp(-io; r ot,2^zí) (8) 

transforms p into the density matrix p in the rotating 
frame as 

p=UpUl (9) 

Note that we can choose the rotation angular velocities 
'j-'roM [i — 1,2) arbitrarily. The time evolution of the 
system is now governed by 

i d £ = [H,pl (10) 
with the transformcd Hamiltonian 

H = UHTJÏ -iU—UÏ = H + H xtl + H I i2- (H) 
dt 

Here the transformed system Hamiltonian is 

dt 

= -(w 0l i - tü Tot< i)I z ® I 

-(CJ 0: 2 - ^Tot,2)I <8>I Z + JIz ® h 

/0 0\ 

J e íAawt 

+ l e - iAu ^ ' 1 ' 
\ 0/ 



where Aw ro t = Wrot,2 — ^rot,i- The transformed control 
Hamiltonians H r fi will be given later. If we take the 
frame co-moving with each spin, which has the angular 
velocities w rot a — üJq,íí the first two terms in Eq. 112|) 
vanish. In the case of heteronuclcus molècules, the con- 
dition |Acl>o| ^> J is always satisfied and thus the ma- 
trix elements in the last line also vanish after averaging 
over time. For example, |Acjo|/27t ~ 400 MHz while 
J/2tt - 200 Hz for 13 C-labeled chloroform at 11 T, for 
which \Au>q\/J ~ 10 6 . Therefore Hq is well approxi- 
mated by 

H a = JI z ®I z . (13) 

When the resonance and co-rotating conditions u) T tj — 
WOjj = w ro t,j are satisfied, the control Hamiltonians in the 
rotating frame 

6tí,í = UH IÍ;i UÍ (14) 
are approximately given as 

fírf.l = —wi t i(co8<fiiI x ®I + sm<piI y ®I), (15) 

fírf,2 = — W1,2(COS02 / ® Ix + SÍn(j)2 I ® Iy) (16) 

after dropping terms rapidly oscillating with freqüències 
2lüoa and Alüq- Note that the factor 2 in front of u>ia 
in Eqs. (@J and © has disappeared in Eqs. (jT5|) and 
(|16fl . This is physically understood as discussed in [ÏÜj : 
a lincarly polarized rf magnètic field oscillating with fre- 
queney cj r f is a superposition of two circularly polarized 
fields with freqüències ±w r f and the effect of the compo- 
nent with — w r f is averaged to vanish. It is also important 
to notice that a pulse with frequeney w T i .i influences only 
the spin i and does not affect the other spin in the rotat- 
ing frame. This is because |Awo| is much larger than the 
inverse of the typical pulse width ~ 1/(10 ps) ~ 100 kHz 
and hence the rf pulse resonating with one spin does not 
have spectral component which affects the other spin. 

In conclusion, the Hamiltonian for a heteronuclcus 
molecule in resonant magnètic fields is 

H = JI Z ® I z 

— wi,i(cos <j>\ I x ® I + sin 0i I y <g) I) 

-UJl3{cOS(j>2l ® Ix + SÍll(j)2 I ® Iy) (17) 

in the rotating frame that has the angular velocities 

B. Homonucleus molècules 

1. Experimental setup 

The NMR setup for homonucleus molècules is shown 
schematically in Fig. [3 Because the difference of the res- 
onance freqüències Alüq — ujq^ — Wo,i is n °t large com- 
pared to luq^ in this case, a common resonance circuit 
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FIG. 2: NMR setup for homonucleus molècules. 



and a power amplifier can be used to control both spins. 
For cytosine in D2O, for example, we find |Awo| /27r ~ 
765.0 Hz while w ,;/27r ~ 500 MHz. Although the differ- 
ence Aojq is small, it still allows us to address respective 
spins individually provided that the pulse width is suffi- 
ciently long. 

The oscillator i generates a continuous rf elèctric wave 
with frequency u) T { t i = Wo,i. The sequencer shapes the 
continuous wave into pulses. When addressing the two 
spins simultaneously, a typical pulse width is of the or- 
der of 10 fxs. On the other hand, when addressing them 
individually, a typical pulse width is of the order of 
27r/ |Acl>o| ~ 1 rns. The rf pulses from the two sequencers 
are mixed and amplified. The coil generates magnètic 
fields and picks up signals from the sample and the re- 
ceiver detects the signals. Due to close resonance freqüèn- 
cies Lüo.i, only one set of resonance circuit and receiver is 
necessary for homonucleus molècules. 



2. Homonucleus molecule in rotating frame 

The Hamiltonian for homonucleus molecule in the lab- 
oratory frame has the identical form to the Hamilto- 
nian for a heteronucleus molecule . Even for homonu- 
cleus molecule the condition |Awq| 3> J is satisfied in 
general. For example, in the case of cytosine in D2O, 
\Alu \/2tt ~ 765.0 Hz while J/2tt - 7.1 Hz, and thus 
the above condition is satisfied. Therefore, the approxi- 
mation used in the derivation of the system Hamiltonian 
(|13|l for a heteronucleus molecule is also applicable to 
derivation of that for a homonucleus molecule. Thus the 
system Hamiltonian of a homonucleus molecule takes the 
form 



H = JI Z ®I Z . 
in the co-rotating frame of each spin. 



(18) 



The control Hamiltonian Hrf i describes the action of 
the resonant magnètic field in the frame rotating with 
angular velocity ui I0 t,i — ^>o,% — ^xï,i- Corresponding 
Hamiltonian H r {^ for homonucleus molecule is consider- 
ably more complicated even when terms rapidly oscillat- 
ing with freqüències 2wo,j and u)o,i + ^0,2 are averaged 
out as 



fírf ,1 = — u\ t \ cos <j)\I x ® I + sin tj>i I y (g) I 
+g cos( Auot + cf>i)I <g) I x 
+g sin( Aw í + <t>\)I <g I y , 



cos 4> 2 I ® I x + sin <p 2 I <g I y ) 



+g 1 cos(-Aw í + 4>2)I X ® I 
+g~ 1 sin(-Aw í + <t>2)Iy ® I 



(19) 



(20) 



If we further assume that the pulse width r are long 
enough so that even slowly oscillating terms in Eqs. I|19|) 
and (|20|l . which contain Ac^o, are averaged out, then 
Eqs. dSJ and lUDI reduce to Eqs. (JT5J and (JTÜJ). Simul- 
taneously, we can tune the pulse width r short enough 
( Jt <C 1) so that the spin-spin interaction l |18| l is negligi- 
ble while pulses arc applicd. Therefore wc conclude that 
an arbitrary pulse sequence designed for heteronucleus 
molècules works for homonucleus ones provided that all 
the hard pulses are replaced by soft pulses whose pulse 
width t satisfies the condition 27t/|Acl>o| < r <C 2n/J. 
We will demonstrate this consequence experimentally in 
the next section, where we set r = 4(27r/ 1 A<Jo|) = 
5.229 ms < 2n/J = 140.8 ms. 



C. Conventional Hamiltonians 

Here we make comparison between the Hamiltonians 
derived in the previous subsection and the Hamiltonian 
for homonucleus molècules used in literature. Conven- 
tionally the system Hamiltonian 



■ff CO nv,0 = -AuJqI (g) I z + JI Z (g> I z 



(21) 



is used to describe spin dynamics in a static magnètic 
field in a rotating frame [15] . Apparently, it differs from 
our Hamiltonians (|18ll . 

We suspect that the Hamiltonian (|21() may be derived 
from the original system Hamiltonian Q via transfor- 
mation from the laboratory frame to the frame rotating 
with a common angular velocity w ro t,i = w ro t.2 = ^0,1 • 
We will show, however, that this choice does not yield 
the Hamiltonian (|21|l in the rotating frame. 

If we take a frame that rotates with a common angu- 
lar velocity equal to wo,i for both spins, transformation 
operator JHJ| becomes 



exp(— iu)o,iIzt) ® Gxp(—iuio,iI z t). 



(22) 
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Then the system Hamiltonian © is transformed into 

d 



H, 



com,0 



U C omHoUç 0m — ÍU C , 
-AüJ I®I z + JI Z 

/ \ 

| 

{ 

\ Ò j 



dt 
I, 



uL 



(23) 



which does not agree with the conventional Hamiltonian 
Another system Hamiltonian in the laboratory framc 
H conv . = -Vo,iI z <8)I- w ,2-í ® I z + JI Z ® I z (24) 

is also sometimes employed in literature 0, 03 but 
this is also different from the original system Hamiltonian 
©. We cannot take the Hamiltonian l|24|) as a correct 
one since we cannot replace ^ fc JI/. £g) Ik by JI Z <g> I z in 
the laboratory frame. 

To illustrate the difFcrcnce between our Hamiltonian 
and the conventional Hamiltonian, let us consider the 
unitary gate 



U E 



exp(— ítïI z ® I z ) 
e ÍTr / 4 



V 












= — ítt/4 



(25) 



which is employed along with one-qubit operations to 
implcment the controlled-NOT gate [ï|. We implement 
the gate Ue from our system Hamiltonian 1181) . as 

Uj(it/J) = cxp(-iirH /J) = cxp(-íttI z <g) I z ). (26) 

In other words, we simply wait for a time interval Tj — 
7r/J without applying any rf pulses. Let us define the 
distance between Ue and Uj(t) as 



\\U E - Uj(t)\\ = y/ti[(U E - U.j(t)y(U E - Uj(t))}. (27) 
This is easily evaluated as 



\U E -Uj(t)\\ = 2v / 2\/l-cos-(Jí-7r). 



(28) 



We observe that the distance vanishes at t = Tj so that 
Uj(Tj) = Ue- We note also that the distance remains 
close to zero in the vicinity t ~ Tj. This robust character 
oïUj (t) was clearly observed in our experiment as shown 
in the next section. 

On the other hand, if we replace Hq in Eq. I|26|) by the 
conventional Hamiltonian l|21(l . the distance between Ue 
and Ujit) becomes 

\\UE-U convJ {t)\\ 
= 2V2^j 1 - cos (^Y^j cos I (JÍ ~ 7r )· (29) 



Therefore, if the conventional Hamiltonian (|21() were 
a correct one to describe the spin dynamics, U conv j(t) 
would not coincide with Ue at í = Tj and the distance 
|| Ue — U conv j(t)\\ should oscillate in the vicinity of the 
t ~ Tj. However, such a ràpid oscillation in time has not 
been observed in our experiment. 



III. EXPERIMENTS 
A. Spectrometer and molècules 

All the data were taken at room temperature with a 
JEOL ECA-500 spectrometer ^8|, where the hydrogen 
Larmor frequeney is approximately 500 MHz. 

We used 0.6 mL, 23 mM sample of cytosine 0| solved 
in D2O. The measured coupling strength is J/2ir = 
7.1 Hz while the frequeney difference is |AíJo|/27r = 
765.0 Hz. The transverse relaxation time T2 is measured 
to be ~ 1 s for both hydrogen nuclei and the longitudinal 
relaxation time T± is ~ 7 s. 

In order to measure the spin states, we apply a reading 
pulse to only one spin, called spin 1, and then obtained 
the spectrum by Fourier transforming the free induction 
decay (FID) signal. The state of the spin 1 is read from 
the sign of the peak in the spectrum while the state of the 
other nucleus (spin 2) is found from the peak position. 



B. Deutsch-Jozsa algorithm 

The Deutsch-Jozsa (DJ) algorithm [2(| is one of the 
simplest quantum algorithms that illustrate the power 
of quantum computation and has been implemented by 
several groups 0, . Let us consider a one-bit func- 
tion / : {0, 1} — > {0, 1}. For a two-qubit register, there 
are only four possibilities for /, whose explicit forms are 
/i(0) = /i(l) = 0, / 2 (0) = / a (l) = 1, / 3 (0) = 0,/ 3 (l) = 
1 and /4(0) = 1,/4(1) = 0. The former two functions 
are said to be "constant" while the latter two are "bal- 
anced" . With the DJ algorithm, we can teli whether a 
given unknown function / is constant or balanced via 
only a single trial. 

Chuang et al. j2j| employed carbon-13 labeled chlo- 
roform, a heteronucleus molecule, as a computational 
resource while Jones and Mosca used cytosine, a 
homonucleus molecule, to execute the DJ algorithm. 
Chuang et al. executed the DJ algorithm using the pulse 
sequences shown in Table ^ According to our previous 
discussions, it should be possible to use the pulse se- 
quences of Chuang et al. for cytosine molècules by simply 
replacing hard pulses with soft ones. 

The results of our quantum computations with cyto- 
sine are summarized in Figs. and We started the 
computation with the thermal equilibrium state since the 
D J algorithm does not require a pure initial state [2ÏJ ■ 
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TABLE I: Control pulse sequences for the Deutsch-Jozsa al- 
gorithm taken from the reference pïj . The functions fi(x) are 
defined in the text. Here, denotes the tt/2 pulses around 
the i-axes (i — x,y,—x,—y). The symbol [n] denotes the n 
pulse around the a;-axis. (1/nJ) denotes the two-qubit entan- 
gling operation produced by turning off the rf pulses during 
the interval 2n/n.J. The pulse sequences are followed by a 
readout tt/2 pulse around the i-axis to the spin 1, which is 
not shown in the Table. The receiver detects the FID (free 
induction decay) signal to measure the spin state. 



Gate 




Pulse sequence 






h 
spin 1 
spin 2 


[|U 


(1/4J) - 
(1/4J) m 


(1/4J) 
(1/47) 


w 


[|]-v 

[|L 


h 

spin 1 
spin 2 


lïlv 

m-v 


(1/4J) - 
(1/4J) m 


(1/4J) 
(1/4J) 




[|]-« 
[f]« 


h 

spin 1 
spin 2 


[fi. 
m-y 


(1/2J) [|] 
[f]» (V2J) [f] 


-V [^]— a 
— 1/ [^]ao 


[f]« 


[f]-« 
[?]* 


h 

spin 1 
spin 2 


[fi» 


(1/2J) [|] 
(1/2J) [f] 


-W [§]— a 


[f]« 


[f]-« 
[f]« 



We varied the J-coupling time interval í in the range 
from 69.8 ms to 71.0 ms. In other words, Au> t/2 was 
swept between 26.7 x 2ir and 27.2 x 2tt. The exact dura- 
tion to produce the designed unitary operator correctly 
is 70.4 ms. We observe from Fig.0th.at the spectra are 
not sensitive to variation of the time interval. Therefore 
we concluded that our Hamiltonian (|18fl accounts for the 
experimental results consistently. 



2. Rf pulse width 

In literature 0, it is recommended to use soft 
pulses whose width is an integral múltiple of 27r/ |Aa>o| 
in order to avoid undesirable effect caused by the term 
AlüqI <E> I z in the conventional Hamiltonian l|21|) . Our 
discussion and experiment show that this tuning is not 
necessary since the relevant Hamiltonian (|18|l does not 
contain the term AlüqI ® I z . 

We executed the DJ algorithms shown in TableQ] with 
different pulse width (5.229 and 6.217 ms). In our set- 
ting, the pulse width 5.229 ms is equal to 4 x 2ir/ \Au>q\, 



1. J-coupling time 

The DJ algorithm employs the J-coupling unitary op- 
erator Uj{t) with t = 2n/4J = Tj/2 or t = 2tt/2J = Tj 
to entangle two spins. The time durations for the two- 
qubit operations are 2tt/4J + 2tï/AJ = Tj for f\ and f% 
and Tj for and f±. Thus the total execution time is 
Tj for all four cases. 

As we discussed when we derived Eq. (|28H , our Hamil- 
tonian i|18fl predicts that Uj(t) does not deviate much 
from the desired unitary transformation even when the 
gate operation time t deviates from the correct value. On 
the other hand, Eq. Q29|l telis us that the conventional 
Hamiltonian l|21l) predicts that Uj(t) sharply depends on 
the timing t and oscillates as cos(Aojoí/2). 

In experiment we executed the DJ algorithm with var- 
ious gate operation time t in the vicinity of Tj and ob- 
served how the resulting spectra depend on t. We em- 
ployed the pulse sequences shown in Table ^ in which 
all hard pulses used by Chuang et al. |2l| were replaced 
with Gaussian soft pulses with the pulse width 5.229 ms. 

The initial state of the molècules is a thermal mixture 
of four states |00), |01), |10), and |11). The DJ algorithm 
does not work when the second qubit is |1) and fails to 
distinguish constant from balanced. On the other hand, 
it works regardless of the state of the first qubit. In Fig. El 
the peaks with a smaller frequency shift (the right peaks) 
are outputs from the initial states 1 01) or |11). In this 
case the algorithm fails to distinguish if /, is constant or 
balanced. The peaks with a larger frequency shift (the 
left peaks) in Fig. are outputs from the initial states 
1 00) or | TO) . In this case the DJ algorithm successfully 
telis us if fi is constant or balanced by the sign of the 
peak (positive for f\ and fi while negative for fs and f±). 




FIG. 3: The FID spectra of the spin 1 in cytosine showing 
output of the DJ algorithms. The sign of each peak indicates 
the state of the spin 1. The location of each peak indicates 
the state of the spin 2. When the initial state of the spin 2 is 
0), it causes a larger shift to the resonance frequency of the 
spin 1 and contributes to a left peak in each curve. Then the 
sign of the left peak discriminates whether fi is constant or 
balanced. The left peak is positive for /i and /2 (constant) 
and negative for fa and fi (balanced). The numbers in the 
left side are durations of the two-qubit operations in ms. The 
correct duration of the two-qubit operations is Tj = 70.4 ms. 
Note that the spectra are insensitive to variation of the two- 
qubit operation time. 
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while 6.217 ms is 4.76 x 27r/ |Aw |. The measured FID 
spectra of the spin 1 are shown in Fig. 0] No significant 
changes in the spectra appeared even if we tuned the 
pulse width to a fractional múltiple of 2ir/ \Aujq\. This 
result proves that the pulse width need not be an inte- 
gral múltiple of 2ir / \Aujq\ to implement a given gate and 
obtain a reasonable spectrum. 



IV. FIELD INHOMOGENEITY 
COMPENSATION 

Here we discuss an experiment to reveal the nature 
of the Hamiltonians 119fl and J3ÜJ, which depict the ac- 
tion of the oscillating magnètic fields on the spins. It is 
common to employ the compensating pulse method to 
suppress errors induced by field inhomogcneity. We will 
show, by employing our Hamiltonian, that the entangling 
operation with the J-coupling is fragile in the prcscncc 
of the compensating pulses and a fine tuning of the gate 
operation time is required. 



A. 7r-pulse pair in J-coupling time 

We have shown in the previous section that it is not 
necessary to tune the J-coupling time very accurately 
since this operation is robust against small change of the 
gate operation time. However, if the system is under the 
influcnce of field inhomogeneity, it may cause an error 
during the J-coupling time. 




FIG. 4: The effect of variation of the pulse widths on the 
spectrum of the DJ algorithms. The pulse widths for spins 
1 and 2 are (a) 5.229 ms and 5.229 ms, (b) 6.217 ms and 
6.217 ms, and (c) 5.229 ms and 6.217 ms, respectively. Ob- 
serve that the spectra are insensitive to the variation of pulse 
width. 



It is well known that this undesired effect caused by 
field inhomogeneity can be compensated by a series of 
hard 7T-pulse pair, whose width is of the order of 10 [xs. 
The best known example may be the CPMG (Carr- 
Purcell-Meiboom-Gill) pulse sequence |l3|. Let us ap- 
ply this technique to Uj(t). Then the pulse sequence for 
Uj(t) is replaced with 

U.u(t) : Uj(t/2) - [tt] - [//(í/2) - [tt], (30) 

where time flows from left to right and [tt] denotes a 
hard pulse that rotates both spins by tt radian. We take 
the rotation axis to be the x-axis of one of the spins, 
for example the spin 1, while the rotation axis for the 
other spin, the spin 2, depends on the time when the 
[tt] is applied, according to the Hamiltonian (|19fl . The 
pulse sequence is represented as the product of unitary 
matriecs 



Uj.(t) = U[ 2 l Uj(t/2) üg Uj(t/2), (31) 



where 

Tj(í) 



-sin(Aw -)/« 



u {2) 



(32) 
(33) 

I x + sm(Aaj t)I <g> Iy}) . 

Note that we put the ratio of resonance freqüències 
9 — w o, 2/^0,1 — 1 f° r the homonucleus molecule. The 
resulting operator C/j*(í) does not coincides with Uj(t). 
The distance between Uj*(t) and C/e is evaluated as 



= exp(— ÍTxI x <X) /) 

x exp(— Í7t[cos(Aüjo — )/ < 

= exp( — íttI x (g) /) 

x exp (— Í7r[cos(Acjot)/ í 



\U E -Uj*(t)\ 
= 2V2W1 



cos — (Jt — tt). 
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(34) 



The distance does not vanish generally because the two 
conditions, cos(AwçJ/2) = ±1 and cos(Jí/4 — 7r/4) = 
±1, are rarely satisfied simultaneously. Moreover, the 
distance is very sensitive to t. 



B. Pseudo-pure state preparation 

An NMR quantum computer must be initialized to the 
pseudo-pure state |00) before executing a specific algo- 



TABLE II: Control pulse sequences to create the pseudo- 
pure state [00) 122^1 . Here, denotes the ^ pulse around 
the i-axes (i = x,y,—x,—y) while (1/2J) denotes the two- 
qubit entangling operation implemented by turning off the rf 
pulses during the interval 2tt/2J. The symbol FG denotes 
application of a pulsed field gradient for spatial labeling. The 
pulse sequence is followed by a readout tt/2- pulse around the 
rz;-axis to the spin 1, which is not shown in the Table. 





Pulse 


sequence 




spin 1 


FG [fj, 


(1/2J) [f] 


_„ FG 


spin 2 


FG 


(1/2J) 


FG 
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However, the spectra were fragile when the duration t 
deviates from the correct value. For example, when the 
duration was set at t = 69.7 ms, the intense signal in- 
dicated that the most of molècules are in the |01) state, 
and not in the desired state 1 00) . 

Thus we conclude that our Hamiltonian i|19|) and i|20l) 
accurately account for the experimental results. 



FIG. 5: The effect of variation of the J-coupling time on 
the spectrum of the pseudo-pure state creation. The spectra 
in the left panel are measured without the 7r-pulse pair in 
the entangling operation. The spectra in the right panel are 
measured with the 7r-pulse pair. The numbers are the two- 
qubit operation time in ms. Note that the spectra generated 
with the 7r-pulse pair are sensitive to variation of the two- 
qubit operation time. 

rithm. This initialization procedure is implemented with 
the pulse sequence [2^ shown in Tablclïïl In this section 
we examine the effect of the compcnsating 7r pulses on 
the initialization procedure. 

The initialization process contains the J-coupling time, 
in which the two spins are entangled by the two-qubit 
operation Uj(t). We vary the gate operation time t to 
introduce an operational error on purpose. It is possible, 
however, to apply a pair of 7r-pulses to compensate this 
error while the J-coupling is under action as instructed in 
Eq. l|3Ufl ■ We have swept the gate operation time between 
69.3 ms and 70.6 ms and the results are summarized in 
Fig.0 

The spectra in the left panel of Fig. [5] were measured 
without applying the 7r-pulse pair in the entangling oper- 
ation. We observed that the spectra are robust against 
small variations of the gate operation time. The in- 
tense peaks with a larger frequeney shift are signals from 
molècules in the |00) state while the smaller peaks are 
error signals from small amount of molècules not in the 
|00) state. 

The spectra in the right panel were measured with the 
7r-pulse pair applied during the entangling operation. We 
observed that the spectra are very sensitive to the vari- 
ation of the J-coupling operation time, although it is 
still possible to adjust the operation time so that the 
desired pseudo-pure state is produced with a good preci- 
sion. When the gate operation time was set at the correct 
value t = 70.3 ms, the spectrum exhibited a sharper peak 
than those in the left panel. This result implies that the 
7r-pulse pair improved the quality of the initialized state. 



V. CONCLUSIONS AND DISCUSSION 

We have derived the relevant Hamiltonian for homonu- 
cleus molècules in NMR quantum computing and shown 
that any pulse sequence for a hctcronucleus molccule may 
be translated into that for a homonucleus molecule by 
simply replacing hard pulses by soft pulses with a prop- 
erly chosen pulse width. We have demonstrated that 
the NMR spectra in several experiments are accounted 
for with our Hamiltonian but not with the conventional 
Hamiltonian found in literature. It was shown in our ex- 
periments that the spectra are robust under small varia- 
tions of the J-coupling operation time as well as of the rf 
pulse widths. Moreover, we provided the theoretical basis 
for field inhomogeneity compensation by a pair of hard 
7r pulses during the entangling operation and verihed it 
experimentally. 

Generalization of the present work to molècules with 
more spins is straightforward. It is easy to find proper 
pulse sequence, either numerically |llj| or by Cartan de- 
composition \WL l23j , once an exact form of the Hamilto- 
nian is obtained. Theoretical analysis as well as exper- 
iments on these subjeets are under progress and will be 
published elsewhere. 
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